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ABSTRACT 

We assess evolution in the Mbh — cr* relationship for quasars in the Sloan Digital Sky Survey- 
Data Release 7 for the redshift range 0.1 < z < 1.2. We estimate the black hole mass, Mbh) using 
the "photoionization method," with the broad H/3 or Mg II emission line and the quasar continuum 
luminosity. For the stellar velocity dispersion, wc use the narrow [O in] or [O n] emission line 
as a surrogate. This study is a follow-up to an earlier study in which we investigated evolution 
in the Mbh — &* relationship in quasars from Data Release 3. The greatly increased number of 
quasars in our new sample has allowed us to break our lowcr-rcdshift subsample into black hole 
mass bins and probe the A/bh — f* relationship for constant black hole mass. The Mbh — < J * 
relationship for the highest-mass (M BH > 10 9 - Mq) and lowest-mass (Mbh < 10 75 M@) black holes 
appears to evolve significantly, however most or all of this apparent evolution can be accounted for by 
various observational biases due to intrinsic scatter in the relationship and to uncertainties in observed 
quantities. The Mbh — c* relationship for black holes in the middle mass range (10 75 < Mbh < 10 90 
Mq) shows minimal change with redshift. The overall results suggest a limit of ±0.2 dex on any 
evolution in the Mbh — c* relationship for quasars out to z m 1 compared with the relationship 
observed in the local universe. Intrinsic scatter may also provide a plausible way to reconcile the wide 
range of results of several different studies of the black hole - galaxy relationships. 
Subject headings: galaxies: active — galaxies: evolution — quasars: general — quasars: emission lines 
— black hole physics 



1. INTRODUCTION 

Strong correlations between the masses of central su- 
permassive black holes and properties of host galax- 
ies, coupled with the seeming ubiquity of such black 
holes, suggest an evolutionary link between black holes 
and their host galaxies. The two strongest correlations 
involve the mass of the black hole, Mbh, and cither 
the host galaxy velocity dispersion, <7* (Gcbhardt et al. 
2000a; Ferrarese & Merritt 2000), or the host galaxy lu- 
minosity (Magorrian et al. 1998) and therefore stellar 
mass. In an effort to determine whether the relation- 
ships hold for all masses and at all times, many studies 
have investigated these relationships for different look- 
back times, covering a wide range of black hole and 
galaxy properties. Over the last decade, the findings 
of these studies have spanned the range of possible re- 
sults. Many works cite positive evolution, in the sense 
that black holes initially appear to outgrow their host 
galaxies (e.g. Peng et al. 2006; Shields et al. 2006; 
Salviander et al. 2007; Woo et al. 2008; Merloni et al. 
2010; Decarli et al. 2010; Targett et al. 2012). Others 
find that black holes are comparatively smaller at high 
redshifts compared with their counterparts in the local 
universe (e.g. Alexander et al. 2008), or that there is no 
significant evolution in the relationship (e.g. Shields et 
al. 2003; Cisternas et al. 2011; Decarli et al. 2012). Var- 
ious scenarios have been suggested for the physical origin 
of the black hole - galaxy relationships, some of which in- 
clude strong evolution (see Portinari et al. 2012 and ref- 
erences therein). A common scenario involves processes 
in which feedback from an active black hole couples with 
star-forming gas in the host galaxy (Silk & Recs 1998; 



Fabian 1999; King 2003; Di Matteo et al. 2005; Mur- 
ray et al. 2005). However, it has been suggested that 
black hole - galaxy relationships have no physical origin 
at all and can arise from a random process of mergers 
of galaxies with uncorrelatcd black hole mass and galaxy 
properties (Jahnke & Maccio 2011). 

Observational bias, whether arising from intrinsic scat- 
ter or from uncertainties in observed properties, present 
a significant difficulty when attempting to assess evo- 
lution in any black hole - galaxy relationship. Several 
observational biases have been identified and modeled 
(e.g. Salviander et al. 2007; Lauer et al. 2007; Shen 
k Kelly 2010; Schulze & Wisotzki 2011), many of which 
can mimic positive or negative evolution in black hole - 
galaxy relationships where no intrinsic evolution exists. 
Volonteri & Stark (2011) note that scatter and observa- 
tional bias can also work to mask an intrinsic negative 
evolution. Lauer et al. (2007) warn that a false signal 
of evolution can be produced when comparing samples 
that have been selected on differing properties. For in- 
stance, local quiescent galaxy samples, which form the 
basis for the black hole - galaxy correlations, are selected 
on galaxy properties, while evolutionary studies are of- 
ten based on samples of active galactic nuclei (AGN), 
which arc primarily selected on properties of the black 
hole. Comparison of black hole and galaxy properties be- 
tween the two samples may yield misleading results. One 
way to mitigate this sort of bias is to compare lower- and 
higher-rcdshift samples selected on the same properties. 

This paper is a follow-up to a previous paper (Sal- 
viander et al. 2007) in which we used quasars from the 
Sloan Digital Sky Survey Data Release 3 (SDSS DR3) 
to assess evolution in the Mbh — c* relationship in the 
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rcdshift range 0.1 < z < 1.2. In this paper, we expand 
our sample by including quasars from the more recent 
SDSS DR7. The greater number of quasars allows us to 
attempt to mitigate some of the effects of biases, par- 
ticularly the comparison bias identified by Lauer et al., 
by assessing the evolution in Mbh — c* over constant 
black hole mass. The Mbh — c* relationship is an im- 
perfect assessor of black hole - galaxy evolution, as it 
appears to be a projection of a more fundamental black 
hole - galaxy relationship involving the effective radius 
of the galaxy bulge (Hopkins et al. 2007); however, it 
is an accessible assessor, particularly for AGN in which 
surrogates from the AGN spectrum can be used in place 
of cr*. 

Giiltekin et al. (2009) give the fit to the local Mbh — o~* 
relationship as 

Mbh = (10 8 - 12 M )(a e /2OO) 4 ' 24 , (1) 

where cr e is the effective velocity dispersion. However, for 
this paper, we adopt the earlier fit given by Tremaine et 
al. (2002) to facilitate comparison with our results from 
Salviander et al. (2007) (hereafter "S07"). Tremaine et 
al. give the fit to the local relationship as 

M BH = (10 8 - 13 M )(a e /2OO) 4 ' 02 , (2) 

All values of luminosity used in this study are 
calculated using the cosmological parameters Hq = 
70 km s" 1 Mpc" 1 , fi M = 0.3, and n A = 0.7. 

2. METHOD 

2.1. Black Hole Masses 

Our method for calculating black hole mass is de- 
scribed in S07 and Shields et al. (2003). Briefly, we 
use the "photoionization method," which assumes that 
the broad line region (BLR) gas orbiting the black hole 
is virialized, such that A/bh = fv 2 R/G. The factor / 
represents a correction to the velocity field for an as- 
sumed geometry for the BLR. The radius of the BLR is 
derived from the radius-luminosity relationship, R cx L 7 , 
calibrated by echo-mapping studies (Wandcl et al. 1999; 
Kaspi et al. 2000, 2005). We adopt the formalism of 
Shields et al. (2003), who chose 7 = 0.5, a choice moti- 
vated by photoionization physics and supported by the 
empirically determined radius-luminosity relationship of 
Bentz et al. (2006). The black hole mass is 

M BH = (10 7 - 69 M Q )vl om L<H, (3) 

where the BLR velocity W3000 is derived from the FWHM 
of the broad H/3 or Mg 11 emission line in units of 3000 km 
s _1 and L44 is the 5100 A continuum luminosity in units 
of 10 44 erg s _1 . For higher redshifts we use the 4000 A 
continuum luminosity calibrated to a 5100 A luminosity 
assuming a power-law function fitted by Vanden Berk et 
al. (2001) for SDSS quasar composite spectra, F v cx v a " 
with a v = —0.44. We assume / = v3/2 for a flattened 
BLR geometry, which is incorporated into equation [3] 

2.2. Stellar Velocity Dispersion 

The stellar velocity dispersion, cr*, is difficult to mea- 
sure directly from the spectra of quasars, so we use emis- 
sion line surrogates from the narrow line region (NLR) 



of the quasar. If the NLR gas orbits in the gravitational 
potential of the host galaxy bulge, this will produce a 
one-to-one proportionality between cr* and the width of 
narrow emission lines. Greene & Ho (2005) examine cor- 
relations between cr* and the widths of [O ill], [O 11], 
and [S II] for a sample of narrow-line AGN, and find 
that tT[o in] tends to be wider than cr* . This is likely due 
to the proximity of the [O ill] -emitting gas to the active 
nucleus producing emission in the blue wing of [O III]. 
Nelson & Whittle (1996) also note that blue wings tend 
to be present on [O ill] emission lines. However, work 
done by Nelson (2000) and Bonning et al. (2005) show 
that in the mean the width of the A5007 [O ill] emission 
line tracks cr* for a range of AGN luminosities. Greene 
& Ho find that the lowcr-ionization lines [O 11] and [S 11] 
track cr* in the mean, though with considerable scatter. 
The correlation between [O in] and [O 11] width shown 
in Figure 4 of S07 further suggests that [O 11] may be 
used as a surrogate for cr* with some confidence. We de- 
fine the surrogate cr* as ctnl = FWHM([0 iii])/2.35 or 
FWHM([0 n])/2.35 for a Gaussian profile. 

3. SAMPLE SELECTION AND SPECTRUM 
MEASUREMENTS 

The quasars for this study were taken from the SDSS 
DR7 (Abazajian et al. 2009). Sample selection and mea- 
surements were carried out as described in S07, except 
as noted below. Briefly, we selected all objects from DR7 
that were classified by the Spectroscopic Query Form as 
quasars. We restricted our search to the redshift range 
0.1 < z < 1.4; the lower limit minimizes the inclusion of 
Seyfert galaxies, while the upper limit is the maximum 
rcdshift at which the [O 11] emission line cleanly appears 
in the SDSS spectral window. 

We carried out spectrum corrections and measure- 
ments using an automated fitting algorithm that first 
corrects the spectra for galactic extinction and rebins 
the wavelength scale from logarithmic to linear with 1.41 
A pixel" 1 . We subtracted off the optical and UV Fe 11 
using templates from Marziani et al. (2003) for the op- 
tical regime and a reconstructed template from Vester- 
gaard & Wilkes (2001) augmented by theoretical data 
from Sigut & Pradhan (2003) for the UV regime. The 
line flux, velocity dispersion (cr GH ), and FWHM of H/3, 
[O in], [O 11], and Mg 11 are measured both before and 
after subtraction of Fe 11. The emission lines are mod- 
eled using Gauss-Hermite functions in order to account 
for any asymmetries (characterized by the /13 parameter) 
and deviations from Gaussianity (characterized by the /14 
parameter). A special case is the FWHM of the [O 11] 
AA3726,3729 doublet, which tends to be unresolved in 
SDSS spectra. We follow the procedure of S07, in which 
we found the [O 11] doublet to be successfully modeled as 
a single line from which the intrinsic width is inferred us- 
ing a calibration curve. Continuum fluxes are measured 
at 4000 A and 5100 A in the observed frame. The instru- 
mental width is subtracted in quadrature from the mea- 
sured line widths, and all measurements are corrected to 
the rest frame. 

We created two subsamples as described in S07: a 
lower-rcdshift sample using H/3 and [O ill] in the red- 
shift range 0.1 < z < 0.81 (hereafter the "H03" sample) 
and a higher-redshift sample using Mg 11 and [O 11] in 
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TABLE 1 

Quality Cut Tolerances for Each Emission Line 



Quality cut H/3 [O III] Mg II [O II] 



EW error 


7% 


5% 


5% 


11% 


FWHM error 


15% 


15% 


15% 


30% 


profile shape (/14) 


0.2 


0.2 


0.2 


0.2 


Reduced x 2 


4 


■1 


4 


4 


BL minimum width (km s" 1 ) 


1500 




1500 





TABLE 2 

Quality Cut Reductions for the 
subsamples 



Quality cut 


# of objects remaining 




H03 


M02 


Initial sample 


22,985 


13,213 


EW error 


14,674 


2214 


FWHM error 


12,270 


2065 


Profile shape (/14) 


11,116 


1334 


Reduced x 2 


10,899 


1322 


BL width minimum 


10,240 


1303 


Visual inspection 


5355 


808 



the redshift range 0.4 < z < 1.35 (hereafter the "M02" 
sample). The H03 and M02 samples consist of 5355 and 
808 objects, respectively. 

We executed a series of quality cuts to remove objects 
with substandard spectra as described in S07 with the 
exception of relaxed standards for the EW and FWHM 
errors when it was found through visual inspection that 
our prior numerical cuts eliminated a number of usable 
spectra. Briefly, we eliminated objects on the basis of 
EW and FWHM measurement errors, excessive deviation 
from Gaussianity (the /14 parameter), reduced x 2 j and a 
minimum width for the broad lines. Table [I] shows the 
quality criteria implemented for each emission line. We 
determined the new EW and FWHM error cuts by select- 
ing a random subsample of 100 DR7 spectra, eliminating 
poor quality spectra through visual inspection alone, and 
adjusting the numerical cuts until we achieved a very 
good match with the hand-selected subsample. There 
are 25,009 DR7 quasars with 0.1 < z < 0.81 and 38,811 
quasars in the range 0.4 < z < 1.4. Table [5] shows the 
number of objects remaining in the subsamplcs after each 
quality cut was implemented. The greatest reduction in 
the number of objects for the M02 sample results from 
the absence of the [O II] emission line from the spectra. 
A great reduction for both subsamples occurs with the 
EW % error cut, which tends to eliminate objects with 
the poorest signal-to-noise ratios. Most of the objects 
eliminated by the EW % error cut were based on poor 
S/N for H/3 or [O 11]. 

We found in S07 that radio loudness did not signifi- 
cantly affect our results, so we have made no effort to 
exclude radio-loud objects from this sample. 

4. RESULTS 

4.1. The Mbh — c* Relationship 

Figure [T] shows the Mbh — 0* relationship for our 
sample of quasars for both the H03 and M02 samples. 



The large red and cyan circles show the mean Mbh and 
°nl (ctio in] or a [o 11] ! respectively) for redshift bins, dis- 
cussed in the following section. The solid line represents 
the local Mbh — &* relationship given by equation [5] and 
is not a fit to the data. The dispersion with respect to the 
local relationship shown in Figure Q] is 0.63 dex and 0.72 
dex for the H03 and M02 samples, respectively. These 
are consistent with the dispersions for the subsamples in 
S07. 

The scatter is mostly attributable to the dispersion in 
the (Tnl- "* surrogacy and the uncertainty in the virial 
estimate for Mbh (see $5]) . Even with the high degree of 
scatter evident in Figure [TJ the mean Mbh and onl do 
center on the local Mbh — cr* relationship for z < 0.45. 
For greater redshifts, the narrow-line width saturates at 
log ctnl ~ 2.24, and the means become increasingly dis- 
placed above the relationship. 
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log <7 NL (km s ») 

Fig. 1. — TheA/BH — o"* relation for the H03 (small open circles) 
and M02 (small closed circles) samples. The large circles show the 
average Mbh ancl °"[0 III] f° r the redshift bins shown in Table [3] — 
red for H03 and cyan for M02. The line is not a fit to data, but 
rather the local A/bh — relationship as given by Tremaine et 
al. (2002). 



4.2. Evolution with Lookback Time 

We compare the photoionization mass with the 
"narrow-line mass," that is, the black hole mass calcu- 
lated with equation |3] using cr[o hi] or n] in place of 
a* (M CT ). The displacement from the local Mbh — c* re- 
lationship is defined as Alog A/bh = log Mbh — log M a . 
If a quasar follows the local Mbh — relationship, it 
will show perfect agreement between the photoionization 
mass and narrow-line mass, i.e. A log A/bh = 0. If 
A log Mbh > 0, this indicates the black hole is pro- 
portionally more massive than would be expected from 
(Tnl and the local Mbh — c* relationship; likewise, if 
A log Mbh < 0, this indicates the black hole is propor- 
tionally less massive than expected from ctnl and the 
local relationship. The mean A log Mbh for the H03 
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sample is +0.13 dex, which heavily weights the abun- 
dant \ow-z quasars in the sample; the mean A log Mbh 
for M02 is +0.59. Both indicate black holes that are 
proportionally more massive than expected. These arc 
consistent with the subsample means of S07. 



1 1 1 1 I 1 1 1 1 1 1 r 




0.5 1 

redshift 

Fig. 2. — A log A/bh — offset from the local Mbh — o", relation — 
as a function of redshift for the H03 (small open circles) and M02 
(small closed circles) samples. The large circles show the average 
A log Mbh f° r the redshift bins shown in Table [3] — red for H03 
and cyan for M02. For most of the redshift bins, the error bars 
are smaller than the marker. 
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log <7 NL (km s" 1 ) 

Fig. 3. — The A'/bh — a * relation for the H03 sample, with colors 
corresponding to the mass bins shown in Figure[4] The large black 
circles show the average Mbh and cr[o hi] for the mass bins. The 
line is not a fit to data, but rather the local AfaH — f* relationship 
as given by Tremaine et al. (2002). 

shows how the intrinsic scatter in the Mbh — &* relation- 
ship leads to systematic offsets from the local relation- 
ship that scale with black hole mass. Note the apparent 
lack of significant evolution in the Mbh — relationship 
except for the highest and lowest black hole masses. 




The pivotal question is whether the A/bh — c* re- 
lationship evolves with lookback time. In Figure [2] wc 
plot A log Mbh as a function of redshift for the H03 
and M02 samples, including the mean for redshift bins 
incremented by Az = 0.1. Tabic |3] shows various quan- 
tities for the redshift bins for H03 and M02. There is a 
general upward trend in A log A/bh with redshift, with 
an overall increase in the average A log Mbh of +0.39 
dex from z = 0.1 to z = 0.8 for H03 only and +0.67 dex 
from z = 0.1 to z m 1.2 for both H03 and M02. Figure 
[2] is consistent with our result from S07. 

Does this upward trend represent real evolution in the 
Mbh — o* relationship with lookback time or is it par- 
tially or wholly the result of other effects contributing to 
the appearance of evolution? We discuss potential biases 
in Sj5j An immediate test is to separate our sample into 
bins restricted by Mbh and determine whether there is 
evidence for evolution with lookback time within each 
of these mass bins. The limited number of objects in 
our sample for S07 prohibited such a test for the higher 
mass bins, but with the inclusion of quasars from DR7 
we have sufficient objects to separate the current H03 
sample into six mass bins, ranging from log Mbh < 7.0 
M Q to log Mbh > 9.0 M Q , that are incremented by 0.5 
dex M Q . (There are too few objects in our M02 sample 
to break it into comparable mass bins.) Figure [3] shows 
the Mbh — c* relation color-coded according to the mass 
bins. Figure [4] shows the same comparison as Figure [2] 
using average quantities for each of the mass bins, and 
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Fig. 4. — A log Mbh — offset from the local A/bh — cr* relation — 
as a function of redshift for the black hole mass bins. 



It is likely that the overall appearance of an upward 
trend in A log Mbh with redshift in Figure [5] is the result 
of a combination of different selection biases. Wc discuss 
these in the following section. 
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TABLE 3 

Average Quantities for Redshift Bins 



z (Number of Objects) vL„(5100) 
erg s — 1 


"NL 
km s _1 


M BH 

M 


A log M B h 


^Edd/ibol 


FWHMbl 

km s _1 


H03 Sample 


0.15 (685) 


43.74 


2.15 


7.60 


+0.07 


-1.01 


3.54 


0.25 (1153) 


43.98 


2.21 


7.75 


-0.02 


-0.92 


3.55 


0.35 (1171) 


44.24 


2.23 


7.92 


+0.07 


-0.84 


3.58 


0.45 (1030) 


44.41 


2.24 


8.02 


+0.15 


-0.76 


3.58 


0.54 (686) 


44.60 


2.24 


8.15 


+0.25 


-0.69 


3.58 


0.64 (419) 


44.76 


2.24 


8.26 


+0.37 


-0.65 


3.60 


0.74 (201) 


44.98 


2.25 


8.39 


+0.48 


-0.55 


3.60 


M02 Sample 


0.47 (43) 


44.62 


2.15 


8.16 


+0.63 


-0.69 


3.59 


0.55 (117) 


44.78 


2.22 


8.22 


+0.43 


-0.59 


3.57 


0.65 (166) 


44.92 


2.20 


8.27 


+0.55 


-0.50 


3.56 


0.75 (177) 


45.02 


2.22 


8.41 


+0.62 


-0.54 


3.61 


0.85 (129) 


45.10 


2.23 


8.40 


+0.55 


-0.45 


3.59 


0.94 (87) 


45.26 


2.22 


8.56 


+0.74 


-0.44 


3.62 


1.04 (61) 


45.35 


2.25 


8.59 


+0.67 


-0.38 


3.61 


1.16 (14) 


45.29 


2.19 


8.49 


+0.83 


-0.35 


3.59 


1.23 (13) 


45.55 


2.23 


8.59 


+0.75 


-0.19 


3.57 



Note. — All quantities are expressed in logarithmic units except for redshift. unl corresponds 
to <T[o in] for the H03 sample and to <j[o II] f° r the M02 sample; FWHMbl corresponds to the 
H/3 FWHM for the H03 sample and to the Mg II FWHM for the M02 sample. Bins with fewer 
than ten objects were excluded (this affected the M02 sample only). 



5. UNCERTAINTIES AND BIASES 

We now discuss potential selection biases that could 
lead to the appearance of evolution in the Mbh — rela- 
tionship with redshift. At the heart of these selection bi- 
ases are scatter and measurement uncertainty. Giiltekin 
et al. (2009) determined an intrinsic rms scatter in the 
Mbh — <y* relation — that is, scatter not due to uncer- 
tainties in measurements and surrogacies — of 0.31 dex 
for early-type galaxies and 0.44 dex for all galaxy types, 
based on the locally-observed sample of galaxies with 
dynamically-determined black hole masses. For AGN 
samples, neither the magnitude of the intrinsic scatter in 
the Mbh — relationship nor the shape of the scatter 
function is known. The 0.63 and 0.72 dex scatter about 
the Mbh~ c* relation for our H03 and M02 samples, re- 
spectively, includes the intrinsic scatter in the Mbh — cr* 
relation, the ~ 0.4 dex uncertainty for virial estimates of 
black hole mass (Vestergaard & Peterson 2006), and the 
0.13 dex scatter in the [O lll]-er* surrogacy (Bonning et 
al. 2005). The combined scatter and uncertainty leads to 
the systematic offsets for the mass bins shown in Figure 
IH with the magnitude of the offsets proportional to the 
scatter. The scatter and uncertainties factor into several 
biases, as discussed below. 

5.1. Malmquist-like bias 

In S07 we simulated the effect of a Malmquist-like bias 
that arises from correlations between quasar luminosity, 
Mbh, and redshift. Galaxies hosting larger black holes 
tend to have higher quasar luminosities and are thus pref- 
erentially selected from the flux-limited SDSS sample. 
Because there is intrinsic scatter in the Mbh — 0* rela- 
tion, a given galaxy in a flux-limited quasar sample will 
typically host a disproportionately large black hole. As- 
suming that quasar luminosity is correlated with Mbh, 
the correlation between luminosity and redshift means 
this selection effect is increased with redshift, and there- 



fore leads to a positive trend in A log A/bh with redshift. 
The strength of this trend is related to the magnitude of 
the scatter. Since it is not known whether the intrinsic 
scatter (both the magnitude and the shape of the scat- 
ter function) in the Mbh — cr* relationship for quasars is 
the same as that for locally-observed galaxies, the degree 
of bias is also not known. Qualitatively, Figure |4] indi- 
cates that this bias is operating in the lowest mass bin, 
in which lower-luminosity objects drop out of the sample 
with increasing redshift (see also Figure [5]) . In our sim- 
ulations of the Malmquist-like bias, scatters of 0.3 dex 
and 0.5 dex produced a factor of 0.19 and 0.44 dex bias, 
respectively, in Mbh for a sample with an average mass 
of log Mbh ~ 8.0, roughly corresponding to the average 
black hole mass of the H03 sample. This range of bias 
is consistent with the overall rise in A log Mbh of +0.39 
dex for the H03 sample (see Figure [2]). 

5.2. Lauer bias 

Lauer et al. (2007) discuss how the appearance of evo- 
lution can arise from comparison of high-redshift samples 
of AGN with local samples of quiescent galaxies, which 
have differing selection criteria. Bias in the Mbh — c* re- 
lationship manifests at high Mbh where there is a dearth 
of high-cr» galaxies to host the biggest black holes. Their 
Figure 3 shows how intrinsic scatter in the Mbh — 0* re- 
lationship gives rise to a higher density of the most mas- 
sive black holes than for the "native" population hosted 
by galaxies with c* prescribed by the Mbh — c* rela- 
tionship. This leaves galaxies with more modest cr* as 
the only remaining homes for the majority of the most 
massive black holes. Such a bias may account for the 
rise in A log Mbh for the highest-mass bin. Lauer et 
al. estimate, for reasonably large scatter, a factor of 3 
bias in Mbh for the Mbh — cr* relation for the biggest 
black holes, consistent with the rise in the largest mass 
bin in Figure [4] Figure [5] shows that, while the black 
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hole mass is roughly constant with redshift for this bin, 
<T[o in] decreases with redshift. A possible explanation 
for this drop off in cr[o hi] with redshift is evolution in 
the <7* function. Chac (2011) compares <7* functions at 
z = and z = 1, based on galaxy data drawn from SDSS 
DR5 and DR6, and shows not only a dramatic decrease 
in large-<7* galaxies at z = 1 compared with z = 0, but 
an increase in the number density of small-cr* galaxies. 
Combined with a a ^-function that already drops off more 
rapidly than the black hole mass function, this leaves in- 
creasingly modest-cr* galaxies to host most of the biggest 
black holes with redshift. 
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Fig. 5. — Log ctnl as a function of redshift (top panel) and log 
Mbh as function of redshift (bottom panel) for the highest-mass 
(log Mbh> 9-0 Mq) H03 bin. See text for discussion. 



5.3. Shen-Kelly bias 

Shcn & Kelly (2010) explore a bias that arises from 
uncertainty in the calculation of Mbh using the pho- 
toionization method. Uncertainties in the photoioniza- 
tion mass arise from, among other quantities, uncertainty 
in the quasar luminosity. Scatter in the luminosity cou- 
pled with the greater abundance of low-mass black holes 
in the black hole mass function means that more of the 
modest black holes are scattered to higher luminosities 
than the other way around. This results in a greater 
number of black holes observed to have higher masses. 
These black holes with observed masses greater than 
their true masses are nonetheless observed to reside in 
galaxies with true cr* corresponding to the true Mbh, 
since the measurements (or proxies) for c* are not sub- 
ject to the same selection biases. Such a bias would give 
rise to the appearance of evolution in the Mbh — c* re- 
lationship. Shen & Kelly estimate a 0.2 - 0.3 dex bias 
in Mbh for £boi > 10 46 erg s _1 , corresponding to the 
average bolometric luminosity for our highest-mass bin. 
This black hole mass bias accounts for at least half of 
the observed offset in this bin at z ~ 0.8. Since this bias 



operates independently of the Lauer bias, it would repre- 
sent an additional contribution to the upward trend for 
the highest mass bin in Figure [4] 

5.4. Schulze-Wisotzki bias 

Schulze & Wisotzki (2011) investigate several biases 
that could affect studies of the evolution in the black hole 

- galaxy relationship in AGN, including the heretofore 
unexplored active fraction bias. Scatter in the black hole 

- galaxy relationship coupled with an active fraction that 
depends on black hole mass can produce either a positive 
or negative bias in a luminosity-limited sample, depend- 
ing on whether the active fraction increases or decreases 
with Mbh- A negative bias occurs if the active fraction 
is inversely proportional to A/bh- In such a case, smaller 
black holes have a greater probability of being fueled, and 
therefore detected in a sample, than larger black holes, 
producing an AGN sample that is more likely to include 
comparatively smaller black holes for a given host galaxy 
property. Schulze & Wisotzki found the active fraction 
bias in a sample like the SDSS to be strongest at the 
luminosity limit of the survey, thus affecting the smallest 
black holes. An active fraction bias adds to any lumi- 
nosity bias present in the sample, and, if present, would 
contribute to either an upward or downward trend for 
the lowest mass bins in Figure 0] depending on the ac- 
tive fraction as a function of A/bh- The bottom panel of 
Figure [6] shows an increasing trend of A/bh for the lower- 
limit of the 7.0 < log Mbh < 7.5 mass bin as a function 
of redshift. (The lowest mass bin is too limited in its 
redshift range to make this analysis.) As shown in the 
top panel of Figure |H1 the downward trend in A log Mbh 
for this mass bin is due to a greater trend of increasing 
a [o in] as a function of redshift. This may be evidence 
of a negative active fraction bias mitigated by a positive 
luminosity bias. 
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Fig. 6. — Log <tnl as a function of redshift (top panel) and 
log Mbh as function of redshift (bottom panel) for the 7.0 < 
log -Mbh < 7.5 Mq H03 bin. See text for discussion. 
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6. DISCUSSION 

6.1. Evolution in the Mbh — C* relationship for SDSS 

quasars 

Is there evidence of genuine evolution in the Mbh — 
cr, relationship or just the appearance of it? Only 
the highest-mass (Mbh > 10 90 Mq) and lowest- mass 
(Mbh < 10 7 5 Mq) bins in Figure [4] show any change in 
A log Mbh with rcdshift. Without knowing the magni- 
tude of the intrinsic scatter or the shape of the scatter 
function for the quasar A/bh — c* relationship, it is not 
possible to determine with accuracy to what degree all 
of these biases affect our sample. However, the obser- 
vational biases and uncertainties described in $5] can ac- 
count for most or all of the apparent evolution in these 
mass bins. The positive trend for the highest-mass bin is 
consistent with both a Lauer selection bias and a Shen- 
Kelly mass-uncertainty bias, while the negative trend for 
the lowest-mass bin is possibly due to a combination of a 
Schulze-Wisotzki active-fraction bias and a Malmquist- 
like luminosity bias in which lower-luminosity quasars 
(and therefore lower-mass black holes) are increasingly 
unlikely to be included in the sample with increasing red- 
shift. The Mbh — &* relationship for black holes in the 
middle mass range (10 7 5 < M B h < 10 9 M ) shows 
minimal evolution with rcdshift. Inspection of the over- 
all results suggests a limit of ±0.2 dex on any evolution 
in the Mbh — c* relationship for quasars from z = 0.1 to 
z = 0.8 compared with the relationship observed in the 
local universe. 



6.2. Reconciling disparate results 

The various studies of evolution in black hole - galaxy 
relationships (see JQ) have painted a confused picture, 
with the conclusions spanning the range of possible re- 
sults. Can these be reconciled into a single, coherent 
picture of black hole - galaxy evolution? 

Lauer et al. (2007) cautioned that comparing AGN 
samples with the local sample of quiescent galaxies could 
produce a false signal of evolution because of the differ- 
ing selection criteria used. Comparing one AGN sample 
with another AGN sample can also be problematic unless 
the selection criteria for both are the same. One way to 
mitigate this problem is to compare lower- z and higher- z 
samples that were selected using the same criteria, which 
we have done to some extent, with the added feature of 
the mass bins. As seen in Figure |4j the magnitude and 
sign of the offset from the local relationship is a func- 
tion of the black hole mass. This may provide some in- 
sight into the problem of the seemingly disparate results 
in studies of black hole - galaxy evolution. Many of the 
AGN studies of evolution focus on black hole - galaxy re- 
lationships not only for relatively small redshift ranges, 
but for a relatively narrow range of black hole masses. 
It is possible that some of these studies are simply ob- 



serving the offset of samples from the local relationship 
that arises from inherent scatter, as opposed to actual 
evolution. 

7. SUMMARY 

We have investigated the Mbh — relationship for 
SDSS DR7 quasars in the rcdshift range 0.1 < z < 1.2 
to assess whether any evolution has taken place. We 
used the "photoionization method" with the widths of 
the broad H/3 or Mg II emission lines and the quasar con- 
tinuum luminosity to calculate black hole masses, and the 
widths of the narrow [O in] or [O n] lines as surrogates 
for the stellar velocity dispersion. We divided our lower- 
redshift H03 sample into black hole mass bins incre- 
mented by 0.5 dex to assess the change in the Mbh — o~* 
relationship as a function of rcdshift for different black 
hole masses. All but one of the mass bins are offset 
from the local relationship, with the magnitude and di- 
rection of the offset determined by the black hole mass 
and the degree of scatter in the relationship. The offsets 
for the lowest-mass and highest-mass black holes appear 
to evolve with redshift; however, various observational 
biases can account for most or all of the apparent evolu- 
tion. The offsets for the middle three mass bins remain 
approximately constant with redshift and suggest a limit 
of ±0.2 dex on any evolution in the Mbh — o~* rela- 
tionship up to z = 0.8. Finally, it may be possible to 
reconcile some of the seemingly disparate results from 
various studies of evolution in the black hole - galaxy re- 
lationships by accounting for offsets in the relationship 
for different masses due to scatter. 



We thank Mark Bottorff for helpful discussions and 
Alyx Stevens for assistance. G.A.S. gratefully acknowl- 
edges the hospitality of Lick Observatory and the support 
of the Jane and Roland Blumbcrg Centennial Professor- 
ship in Astronomy. 

Funding for the Sloan Digital Sky Survey (SDSS) has 
been provided by the Alfred P. Sloan Foundation, the 
Participating Institutions, the National Aeronautics and 
Space Administration, the National Science Foundation, 
the U.S. Department of Energy, the Japanese Monbuk- 
agakusho, and the Max Planck Society. The SDSS 

The SDSS is man- 



Web site is http://www.sdss.org/ 
aged by the Astrophysical Research Consortium (ARC) 
for the Participating Institutions. The Participating 
Institutions are The University of Chicago, Fcrmilab, 
the Institute for Advanced Study, the Japan Participa- 
tion Group, The Johns Hopkins University, the Korean 
Scientist Group, Los Alamos National Laboratory, the 
Max-Planck- Institute for Astronomy (MPIA), the Max- 
Planck- Institute for Astrophysics (MPA), New Mexico 
State University, University of Pittsburgh, University 
of Portsmouth, Princeton University, the United States 
Naval Observatory, and the University of Washington. 



REFERENCES 



Abazajian, K. N., Adelman-McCarthy, J. K., Agiieros, M. A., et 

al. 2009, ApJS, 182, 543 
Alexander, D. M., Brandt, W. N., Smail, I., et al. 2008, AJ, 135, 

1968 

Bentz, M. C, Peterson, B. M., Poggc, R. W., Vestergaard, M., & 
Onken, C. A. 2006, ApJ, 644, 133 



Bonning, E. W., Shields, G. A., Salviandcr, S., & McLure, R. J. 

2005, ApJ, 626, 89 
Chae, K.-H. 2011, MNRAS, 413, 887 

Cisternas, M., Jahnke, K., Bongiorno, A., et al. 2011, ApJ, 741, 
Lll 

Dccarli, R., Falomo, R., Treves, A., et al. 2010, MNRAS, 402, 
2453 



8 



Salviandcr & Shields 



Dccarli, R., Falomo, R., Kotilainen, J. K., et al. 2012, Advances 

in Astronomy, 2012 
Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 

604 

Fabian, A. C. 1999, MNRAS, 308, L39 
Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 
Gcbhardt, K., Bender, R., Bower, G., et al. 2000, ApJ, 539, L13 
Greene, J. E., & Ho, L. C. 2005, ApJ, 627, 721 
Giiltckin, K., Richstone, D. O., Gcbhardt, K., et al. 2009, ApJ, 
698, 198 

Jahnke, K., & Maccio, A. V. 2011, ApJ, 734, 92 

Kaspi, S., Smith, P. S., Netzer, H., Maoz, D., Jannuzi, B. T., & 

Giveon, U. 2000, ApJ, 533, 631 
Kaspi, S., Maoz, D., Netzer, H., Peterson, B. M., Vestergaard, 

M., & Jannuzi, B. T. 2005, ApJ, 629, 61 
King, A. 2003, ApJ, 596, L27 

Lauer, T. R., Tremaine, S., Richstone, D., & Faber, S. M. 2007, 
ApJ, 670, 249 

Magorrian, J., Tremaine, S., Richstone, D., et al. 1998, AJ, 115, 
2285 

Marziani, P., Sulentic, J. W., Zamanov, R., Calvani, M., 
Dultzin-Hacyan, D., Bachcv, R., & Zwitter, T. 2003, ApJS, 
145, 199 

Merloni, A., Bongiorno, A., Bolzonella, M., et al. 2010, ApJ, 708, 
137 

Murray, N., Quataert, E., & Thompson. T. A. 2005, ApJ, 618, 569 

Nelson, C. H. 2000, ApJ, 544, L91 

Nelson, C. H., & Whittle, M. 1996, ApJ, 465, 96 



Peng, C. Y., Impey, C. D., Rix, H.-W., et al. 2006, ApJ, 649, 616 
Portinari, L., Kotilainen, J., Falomo, R., & Decarli, R. 2012, 

MNRAS, 420, 732 
Salviander, S., Shields, G. A., Gcbhardt, K., & Bonning, E. W. 

2007, ApJ. 662, 131 
Schulze, A., & Wisotzki, L. 2011, A&A, 535, A87 
Shen, Y., & Kelly, B. C. 2010, ApJ, 713, 41 

Shields, G. A., Gebhardt, K., Salviander, S., Wills, B. J., Xie, B., 
Brotherton, M. S., Yuan, j., & Dietrich, M. 2003, ApJ, 583, 
124 (S03) 

Shields, G. A., Menezes, K. L., Massart, C. A.. & Vandcn Bout, 

P. 2006, ApJ, 641, 683 
Sigut, T. A. A., & Pradhan, A. K. 2003, ApJS, 145, 15 
Silk, J., & Rees, M. J. 1998, A&A, 331, LI 
Targett, T. A., Dunlop, J. S., & McLure, R. J. 2012, MNRAS, 

420, 3621 

Tremaine, S., Gebhardt, K., Bender, R., et al. 2002, ApJ, 574, 740 
Vanden Berk, D. E., et al. 2001, AJ, 122, 549 
Vestergaard, M., & Peterson, B. M. 2006, ApJ, 641, 689 
Vestergaard, M., & Wilkes, B. J. 2001, ApJS, 134, 1 
Volonteri, M., & Stark, D. P. 2011, MNRAS, 417, 2085 
Wandel, A., Peterson, B. M., & Malkan, M. A. 1999, ApJ, 526, 
579 

Woo, J.-H., Treu, T., Malkan, M. A., & Blandford, R. D. 2008, 
ApJ, 681, 925 



